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Model  predictions  of  bound  (i.e.,  nonlinearly  forced  by  and  coupled  lo  wave  groups)  infragravity  wave 
enei^y  are  compared  with  about  2  years  of  observations  in  8-  to  13-m  depths  ai  Imperial  Beach.  California,  and 
Barbers  Point,  Hawaii.  Frequency^itectional  spectra  of  free  waves  at  sea  and  swell  frequencies,  estimated  with 
a  small  array  of  four  pressure  sensors,  ate  used  to  predict  the  bound  wave  spectra  below  0.04  Hz.  The  predicted 
total  bound  wave  energy  is  always  leas  than  the  observed  infragravity  energy,  and  the  underprediction  increases 
with  increasing  water  depth  and  especially  with  decreasing  swell  energy.  At  most  half,  and  usually  much  less, 
of  the  observed  infragravity  energy  is  bound.  Bound  wave  spectra  are  alto  predicted  with  data  from  a  tingle 
wave  gage  in  183-m  depth  at  Point  Conception,  California,  arid  the  assumption  of  unidirectional  sea  and  swell. 
Even  with  energetic  swell,  less  than  10%  of  the  total  observed  infragravity  energy  in  183-m  depth  is  bound. 
Free  waves,  either  leaky  or  edge  waves,  are  mote  energetic  than  botmd  waves  at  both  the  shallow  and  deep  sites. 
The  low  level  of  infragravity  energy  observed  in  183-m  depth  compared  with  8-  to  I3-m  depths,  with  similarly 
moderate  sea  and  swell  energy,  suggests  that  leaky  (and  very  high-mode  edge)  waves  contnbule  less  than  10% 
of  the  ittfragravity  energy  in  8-13  m.  Most  of  the  free  infragravity  energy  in  shallow  water  is  refraciively 
tra^l^  and  teacli|le^^aj^ 


1.  Introduction 

Infragravity  waves  are  believed  to  be  an  important  factor  in 
several  nearshore  processes.  The  purpose  of  this  paper  is  to  esti¬ 
mate  the  contribution  of  bound  waves  to  infragravity  energy 
observed  well  outside  the  surf  zone  in  depths  of  both  ~10  and 
~200m.  Infragravity  motions  (typical  periods  of  25-200  s  on 
Pacific  coasts)  coupled  to  incident  wave  (typical  periods  of 
4-25  s)  groups  were  first  observed  in  roughly  15  m  depth  by 
Marik  [1949]  and  Tucker  [1950],  who  showed  suggestive  correla¬ 
tions  between  wave  groups  and  low-frequency  motions.  In  both 
cases  the  infragravity  wave  heights  were  about  10%  of  the 
incident  wave  heights. 

Weak  nonlinear  interactions  between  first-order  free  waves  (sea 
and  swell)  of  nearly  equal  frequency  is  one  possible  mechanism  of 
generating  infragravity  waves  bound  (i.e.,  phase  coupled)  to 
groups  of  higher-frequency  waves  [Diesel,  1952;  Longuet-l/iggins 
and  Stewart,  1960,  1962;  llasselmann,  1962;  and  others).  Phase 
coupling  between  infragravity  waves  and  sea  and  swell  has  been 
observed  in  8-  to  18-m  depths  [llasselmann  el  al.,  1963;  Meadows 
el  at.,  1982;  Elgar  and  Guza,  1985;  Elgar  et  al.,  1989],  Biphases 
of  pressure  in  depths  of  8  m  [Elgar  et  al.,  1989],  11m  [llassel¬ 
mann  et  al.,  1963],  and  18  m  [Elgar  and  Guza,  1985]  sliow  that 
interactions  between  waves  near  the  peak  frequency  of  the  spec¬ 
trum  drive  infragravity  motions  with  a  depression  of  sea  level 
under  the  largest  waves  in  a  group,  consistent  with  bound  wave 
theory  for  nearly  colHnear  waves.  These  observations  show  that 
bound  waves  arc  detectable,  but  quantitative  comparisons  of 
predicted  bound  wave  and  observed  infragravity  energy  levels  ore 
surprisingly  limited.  Using  two  70-mm-long  records  of  energetic 
wave  conditions  in  40-m  depth.  Sand  [1982i]  found  very  good 
agreement  (predictions  within  about  5%  of  the  total  observed 
infragravity  energy)  at  a  site  in  the  North  Sea  a  few  hundred 
kilometers  distant  from  the  coast,  with  the  assumption  that  the 
bound  waves  and  forcing  waves  (sea  and  swell)  were  collincar. 
However,  other  observations  suggest  that  free  infragravity  waves, 

well  as  bound  waves,  are  important. 
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Following  the  early  observations  by  Munk  [1949)  and  Tucker 
[1950]  suggesting  that  the  seaward  propagating  infragravity  wave 
amplitude  was  at  least  as  large  as  the  amplitudes  of  the  shoreward 
propagating  bound  wave,  Longuet -Higgins  and  Stewart  (1962) 
speculated  lhat  the  incoming  bound  wave  somehow  reflects  from 
the  shoreline  and  radiates  seaward  as  a  free  wave.  Numerical 
models  [Symonds  et  al.,  1982]  suggest  that  slow  modulation  of  the 
breakpoint  position  at  the  group  frequency  results  i.i  long-wavc 
radiation  seaward  from  the  breakpoint,  but  laboratory  results  are 
inconclusive.  Kostense  [1984]  measured  the  amplitudes  of 
infragravity  waves  in  shallow  water  induced  by  wave  grouping  in 
a  long  wave  channel  with  a  plane  beach  at  one  end.  The  observed 
and  theoretical  bound  wave  amplitudes  agreed,  but  there  were 
significant  differences  between  the  experimentally  observed  out¬ 
going  free  wave  amplitudes  and  those  predicted  by  Symonds  et  al. 
[1982],  Kostense  [1984]  attributed  these  differences  to  assump¬ 
tions  made  in  the  Symonds  et  al.  model  which  failed  to  include  the 
incoming  bound  wave  and  assumed  complete  reflection  at  the 
shoreline,  both  assumptions  affect  the  outgoing  wave  amplitude. 
However,  similar  laboratory  experiments  by  Mansard  and  Bartkel 
[1984]  do  not  appear  to  show  the  presence  of  breakpoint-forced 
outgoing  waves.  It  is  utKlear  how,  or  if.  shoreward  propagating 
bound  wave  energy  is  reflected  from  very  shallow  water. 

Another  generation  mechanism  of  infragravity  motions  has 
been  proposed  by  Gallagher  [1971],  who  showed  theoretically 
that  certain  directional  distributions  of  the  incident  wave  field  can 
resonantly  excite  edge  waves  (free  waves  tra(^d  in  shallow 
water  by  reflection  and  refraction).  Laboratory  results  confirm 
that  these  directional  distributions  indeed  yield  elevated  infragrav- 
ily  energy  levels  [Bowen  and  Guza,  19781.  Recent  observations 
with  arrays  of  flow  meters  in  very  shallow  depths  (<3  m)  show 
that  infragravity  energy  levels  increase  substantially  near  the 
shoreline  and  include  a  signifleant  contribution  from  edge  waves 
[Huntley  et  at..  1981;  OUman-Shay  and  Guza,  1987;  Hawd  et  al., 
1991J.  Roughly  30-50%  of  the  run-up  variance  at  infragravity 
frequencies  is  estimated  to  be  contributed  by  edge  waves  with 
mode  numbers  i  2  [Ollman-Shay  and  Guza,  1987].  While  the 
presence  of  low-mode  edge  waves  cm  explain  the  observed 
energy  in  the  longshore  velocity  fleld,  low-mode  edge  waves 
alone  cannot  describe  the  observed  cross-shore  velocity  field. 
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Other  waves  affecting  the  cross-shore  velocity  much  more  than 
the  longshore,  such  as  high-mode  edge  waves  and/or  leaky  waves, 
are  clearly  present  in  field  data  from  surf  zones. 

Although  the  generation  mechanism  of  free  infragravity  waves 
is  unclear,  bound  waves,  edge  waves,  and  reflections  of  <9  (1)  at 
the  shoreline  (producing  standing  waves)  have  been  detected  in 
many  field  observations  in  very  shallow  water  (<  3  m)  within  a 
few  hundred  meters  of  the  shoreline  [Suhayda,  1974;  Huntley, 
1976;  Sasaki  et  al.,  1976;  Sasaki  and  Horikawa,  1978;  Holman, 
1981;  Huntley  el  al.,  1981;  Guza  et  al.,  1984;  Huntley  and  Kim, 
1984;  Guza  and  Thornton,  1985;  Elgar  and  Guza,  1985:  List, 
1986;  and  others].  Preliminary  results  from  several  hours  of  data 
from  an  array  further  offshore  (8-m  depth,  about  1  km  from  shore) 
suggests  that  high-mode  edge  waves  (mode  number  23)  may 
sometimes  contribute  signific<urtly  to  (or  dominate)  the  infragrav¬ 
ity  wave  field  well  outside  the  surf  zone  where  low-mode  edge 
waves  are  not  energetic  lOltman-Shay  et  al.,  1989;  Elgar  et  al., 
1989;  J.  Oltman-Shay,  personal  communication,  1990],  However. 
Sand's  [19826]  results  (based  on  the  assumption  of  unidirectional 
waves)  suggest  that  bound  waves  dominate  in  40  m  depth.  In 
addition,  directional  spectr®  measured  in  12.5-m  depth,  approxi¬ 
mately  6  km  off  the  New  Jersey  coast,  show  the  infragravity 
waves  (15-  to  85-s  periods)  to  be  propagating  onshore  in  the  same 
direction  as  the  swell  (15-s  period).  Seaward  propagating  leaky 
waves  or  longshore  propagating  edge  waves  were  not  detectable 
[Goodman  et  al.,  1989]. 

The  relative  importance  of  bound  waves,  high-mode  edge 
waves,  and  leaky  waves,  as  a  function  of  position  on  the  continen¬ 
tal  shelf  and  other  variables,  is  not  known.  Here  we  assess  the 
relative  importance  of  bound  waves  well  outside  the  surf  zone  in 
the  depth  range  8-13  m  (0.4-0.8  km  offshore),  where  bound  wave 
theory  has  been  applied  with  qualitative  success  [Hasselmann  et 
al.,  1963].  This  data  set,  which  also  includes  sensors  in  4.5-  and 
7.0-m  depths,  consists  of  2  years  of  observations  from  Barbers 
Point,  Hawaii,  and  Imperial  Beach,  California.  Several  days  of 
data  from  183-m  depth,  10  km  offshore  on  the  continental  shelf 
near  Point  Conception,  California,  are  also  considered. 

Bound  wave  theory  is  reviewed  in  section  2.  Examples  are 
given  which  illustrate  the  strong  dependence  of  bound  wave 
energy  levels  on  directional  spreading,  water  depth,  and  other 
parameters.  The  field  sites  and  data  sets  are  described  in  section 

3.  Bound  wave  model-data  comparisons  are  presented  in  section 

4.  We  conclude  that  bound  waves  contribute  less  than  half  of  the 
observed  infragravity  energy  in  8-  to  13-m  depths,  and  a  much 
smaller  fraction  in  183-m  depth.  The  remainder  is  free  (either 
leaky  or  edge)  wave  energy.  In  section  5  we  show  that  only  a 
small  amount  (at  most  10%)  of  the  free  wave  energy  in  8-  to  13-m 
depth  reaches  the  gage  in  '83-m  depth.  Most  of  the  free 
infragravity  energy  in  8-  to  13-m  depth  is  refractively  trapped 
edge  waves.  The  results  are  summarized  in  section  6. 


2.  Bound  Wave  Theory  Review 

Predictions  of  infragravity  bound  waves  in  constant  depth  are 
obtained  by  expanding  the  inviscid,  incompressible,  irrotational 
wave  equations  and  boundary  conditions,  by  the  method  of 
Stokes,  to  second  ordw  in  the  wave  slope  [Diesel,  1952;  Longuet- 
Higgins  and  Stewart,  1960,  1962,  1964;  Hasselmann,  1962;  and 


others].  At  lowest  order  the  sea  surface  elevation  is  assumed  to  be 
a  linear  sum  of  free  waves  (sea  and  swell), . ' 

la  It  *■  '■* 
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where  a}  =  gk, tanhEi h,  h  is  the  depth,  (x  , y )  are  the  horizontal 
coordinates,  and  Sus,  SB  ate  the  frequency  and  dir^tional  resolu¬ 
tion,  respectively.  A,^  and  are  the  first-order  free  wave 
amplitude  density  and  phase,  respectively,  of  waves  with  radian 
frequency  co,(=nSu)),  vector  wavenumber  t,(k.  =  j).  and 

propagation  direction  0, .  The  lowest  and  highest  free  wave  (sea 
and  swell)  frequencies  are  nusSia  and  nst&u ,  respectively,  and 
=  360/ SB  is  the  number  of  directional  bands  at  each  frequency. 

The  second-order  (bound  wave)  sea  surface  elevation  with  only 
the  difference  frequencies  retained  are 

r?«>(W>)="S  I 

n  =  ^  =  I  r  =  l 

•  cos[Aib;cos0a  +  Akysinfla  -cojt  +  &^]Sa)SB  ,  (2) 

where  the  bouixl  wave  frequency  a>j  =0),,j  -tu,.  the  bound 

wave  vector  wavenumber  aI  =  .y,  -  =  (Akcosda  , 

Aksinfia),  Ak  =  |A)i  |,  aixl  £tp  =0, -^>4  is  the  bound  wave 
phase.  The  bound  wave  sea  surface  elevation  coupling  coefficient 
,C  is  a  function  of  the  free  wave  frequencies  et),,y  ,  ft),  and  the 
angular  difference  A6  =  0,  .  Note  that  the  bound  wave 

propagation  direction  0*  is  determined  from  the  vector  difference 
Sk  and  is  not  equal  to  A0. 

Assuming  that  the  first-order  free  waves  have  random  phases, 
the  infragravity  bound  wave  sea  surface  elevation  frequency  spec¬ 
trum  ,5<b(ft)y)  is  [Hasselmann,  1962] 

.S®(ft)y)  =  2&)  I  I  I 

»  =  fly  V  '  \ 

•S(')(O).,0,>S('>(fi)„y,O,)»J.  (3) 

where  the  free  wave  (sea  and  swell)  frequency-directional  ^xsc- 
trum  is  5*'*(ft),  .0, )  and  the  bound  wave  energy  at  frequency  (Oj 
has  been  summed  over  all  the  infragravity  wave  directions  0s  . 
The  bound  wave  pressure  spectrum  is  given  by  equation  (3)  with 
,C  replaced  by  pC ,  the  bound  wave  pressure  coupling 
coefficienL  Frequency-directional  wavenumber  spectra  of  botmd 
waves  are  obtained  with  separate  summations  over  the  free  wave 
components  contributing  to  specific  infragravity  frequencies  and 
Wavenumbers.  The  bound  wave  sea  surface  elevation  spectrum 
is  related  to  the  smoothed  instantaneous  wave  energy  his¬ 
tory  (SrWEH)  spectrum  [Funke  and  Mansard,  1979;  Sand, 
1982b].  The  SIWEH  spectrum  describes  the  grouping  structure  of 
the  sea  and  swell  wave  field  but  lacks  the  coupling  coefficient  ,  C 
theoretically  relating  the  amplitudes  of  groups  and  bound  waves, 
and  can  thus  only  very  qualitatively  mimic  the  physics  of  bound 
wave  theory. 

The  coupling  coefficients  ,C  arxl  ,C  (at  vertical  distano:  z 
measured  from  the  sea  surface)  are  given  by  [Hasselmann,  1962; 
Sand,  1982a] 
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where  the  sign  of  the  coupling  coefficient  relates  the  phases  of  the 
bound  waves  and  the  free  wave  groups.  When  the  interacting  free 
waves  are  collinear  (dt?  =  0°),  the  magnitude  of  the  sea  surface 
elevation  coupling  coefficioit  |<C)  is  larger  than  the  bottom 
(z  = -A )  pressure  coupling  coefficient  \fC  \  in  all  depths.  How¬ 
ever,  in  water  shallow  for  the  sea  ana  swell.  ia)}hlg  «  1) 
1,C  1  ~  \pC  I,  whereas  in  deep  water,  \,C  I  y>  |,C  |  (Figure 
la).  For  directionally  spread  (60  *  0”)  incident  waves  there  are 
intermediate  depths  and  incident  wave  frequencies  where  |pC  |  is 
slightly  larger  than  |,C|  (Figure  lb,  60  =  20°),  and  the  bound 
wave  pressure  fluctuations  at  the  seafloor  are  actually  larger  than 
ail  unattenuated  hydrostatic  response  to  infragravity  sea  surface 
displacements.  The  bound  wave  response  varies  somewhat  with 
changes  in  the  ratio  of  bound  to  free  wave  frequency  ((Ojla, ;  Fig¬ 
ure  la),  and  strongly  as  a  function  of  the  angular  spread,  depth, 
and  free  wave  frequency  (Figure  lb).  In  shallow  water  (for  the 
sea  and  swell)  the  magnitude  of  the  coupling  coefficients  are  larg¬ 
est  when  the  interacting  waves  are  collinear  (60  =  0®),  and  bound 
waves  foreed  by  even  slightly  directionally  spread  seas  are 
theoretically  significantly  smaller  than  with  unidirectional  seas 
(compare  60-0°  with  60  =  10°  in  Figures  lb  and  2a  and  Table 
1 ).  In  contrast,  when  the  sea  and  swell  waves  are  in  deep  water, 
the  sea  surface  elevation  coupling  coefficient  has  Che  largest  mag¬ 
nitude  for  large  angular  spread  (Figure  lA),  and  free  waves  with  a 
directionally  broad  spectrum  may  force  more  bound  wave  energy 
at  the  sea  surface  than  is  forced  by  a  relatively  narrower  spectrum 
(Figures  2c  and  2d  and  Table  1). 

For  collinear  free  waves,  the  bound  infragravity  wave  is  180° 
out  of  phase  with  the  free  wave  grotips  (i.e.,  ,C  and  pC  are 
negative)  and  the  infragravi^  wave  trough  occurs  beneath  the 
high  waves  in  the  group  ILongu^l-f/iggiru  and  Stewart,  1962]. 
For  moderately  large  60(i3O°  in  deep  water:  Figure  Ih),  the  sea 
surface  elevation  coupling  coefficient  changes  sign  and  the  analo¬ 
gous  phase  is  0°.  With  increasingly  large  60,  both  groups  and 
bound  waves  have  progressively  shorter  spatial  scales  in  both  hor¬ 
izontal  directions.  When  the  angular  spread  between  interacting 
free  waves  of  about  the  same  frequency  is  very  large  (60  >  90°), 
the  bound  wavenumber  is 

\St\  =  \t.,jp-Jt,^\>'i2\t..iA  (5) 

and  the  infragravity  bound  wave  is  shorter  than  the  forcing  free 
waves.  These  bound  waves  are  coupled  to  the  "checkerboard" 
pattern  of  short-crested  waves  rather  than  the  "group"  structure  of 
long-crested  waves.  The  bound  wave  crest  occurs  UTtder  maxima 
in  the  checkerboard. 


It  does  not  seem  to  have  been  previously  noted  that  with  very 
oroad  directional  distributions  in  deep  water,  these  short- 
wavelength  waves  can  theoretically  dominate  the  infragravity 
bound  wave  field  at  the  sea  surface,  owing  to  the  maximum  of 
,C  at  very  large  60.  They  will  not,  however,  penetrate  deeply 
into  the  water  column.  These  short-wavelength,  low-frequency 
waves  are  generated  by  the  same  interacting  waves  (i.e.,  direction- 
ally  opposing  wave  trains  of  nearly  equal  frequencies)  that  gen¬ 
erate  the  long-wavelength,  high-frequency  motions  that  dominate 
the  pressure  spectrum  on  the  seafloor  in  deep  water  and  excite 
microseisms  [Longuet-Higgins,  1950].  In  very  deep  water. 
O  (5000  m),  the  longest  infragravity  bound  waves,  driven  by  col- 
linear  interactions,  are  short  [AA  =(0)}^j  -(oi)/g]  compared  with 
the  water  depth,  attenuation  at  the  seafloor  is  extreme,  and  the  bot¬ 
tom  pa'essure  field  is  dominated  by  free  infragravity  waves 
[wavenumber  (0„j  -(0,'^lg  \  radiated  frmn  coastlines  with  ener¬ 
getic  sea  and  swell  [Webb  et  al.,  1991].  The  most  energetic  free 
mfragravity  waves  observed  in  the  deep  Pacific  [WeW  et  al., 
1991]  have  variance  of  about  0.06  cm^.  However,  at  the  sea  sur¬ 
face  the  predicted  bound  wave  variance  is  -0.2  cm^  for  energetic 
(variance  -l(f*  cm^)  sea  and  swell  (Table  1).  Thus  infragravity 
bound  wave  energy  levels  might  be  comparable  to  (or  larger  than) 
free  infragravity  waves  at  the  sea  surface  in  tlie  deep  ocean,  even 
though  free  waves  dominate  the  bottom  pressure. 

As  nearly  collinear  first-order  free  waves  (sea  and  swell)  enter 
very  shallow  water,  near-resonances  occur  because  the  sea  and 
swell  group  velocity  approaches  the  phase  velocity  of  firee 
infragravity  waves.  In  shallow  water,  ,C^  and  become  very 
large  for  small  60  (Figure  1),  the  assumption  (of  bound  wave 
theory)  that  predicted  bound  waves  have  much  smaller  amplitudes 
than  the  first-order  free  waves  is  violated  (Figure  3a),  and  qualita¬ 
tive  results  at  best  can  be  expected  [Longuet-Higgins  and  Stewart, 
1964].  Note  that  the  depth  at  which  bound  wave  spectral  dettsity 
exceeds  the  first-order  free  wave  densi^  depends  on  both  the 
first-order  wave  amplitudes  and  frequencies.  If  the  swell  energy 
in  Figure  3a  were  reduced  by  a  factor  of  10,  the  predicted  bound 
wave  energy  would  decrease  by  10^  (equation  (3)).  The  shallow 
water  limit  of  applicability  for  bound  wave  theory  is  not  well 
understood. 

Fot  simplicity,  consider  the  sea  surface  elevation  coupling 
coefficient  for  collinear  free  plane  waves,  with  both  the  free  and 
bound  waves  in  shallow  water,  to  0(kh  )^ 

.C(shallow)  =  -i-rTrT —  =  -4-  .  .g - ^  (6) 
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Fig.  1.  The  nondimeiuianal  (note  the  factor  of  h)  bound  wave  sea  surface  elevation  coupling  coefficient  (|,CA  | ,  solid  and  dashed 
lines)  and  bottom  (r  =  -A)  pressure  coupling  coefficient  (|pCA  | ,  dotted  lines)  as  a  functitxt  of  the  first -order  free  wave  radian  fre¬ 
quency  CO,,  bound  wave  radian  frequency  (By.  depth  k  ,  angular  separation  A0,  and  ratio  tuy/o),  (after  Oiltsen  Hansen,  1978; 
Send,  19824].  (n)  Various  values  of  tuyfai,  and  constant  Ad  w  0.  (h)  Various  values  of  Ad  and  constant  (Uy/tu,  =  0  1.  Solid  and 
doited  lines  indicate  negative  values  of  ,  C  and  pC  ,  respectively,  and  dashed  lines  arc  for  positive  values  of  ,  C . 


ami  the  ratio  of  bound  to  free  wave  energy  is  approximately  pro¬ 
portional  to  Ur\  where  Ur  ={alh)Hkh)'^  is  the  Ursell  number 
and  a  and  k  are  characteristic  free  wave  amplitude  and 
wavenumber,  respectively.  Thus  bound  wave  energy  in  shallow 
water  increases  rapidly  with  decreasing  first-order  free  wave  fre¬ 
quencies  and  decreasing  depth.  As  shown  in  Figure  3,  given 
e<]ually  energetic  unidirectional  narrow  banded  swell  (0(0.07 
Hz))  and  sea  (0(0.2  Hz))  in  the  same  (10  m)  depth,  the  swell  will 
force  about  10^  times  more  bound  wave  energy  below  0.04  Hz 
chan  the  sea. 

Nelson  et  cd.  [1988]  compared  infragravtty  wave  heights  (//“*") 
measured  in  12  m  with  sea  and  swell  wave  heights  (//js)  and  fre¬ 
quencies  (,(Oss )  measured  in  SO  m  depth  and  developed  an  empiri¬ 
cal  equation  which  can  be  written  //■*•  -  (//sj /fflsi  I'  ’’ .  a  trend  in 
qualitative  agreement  with  the  shallow  water  limit  of  bound  wave 
theory  (bound  wave  height  -  (llss/Ofss)^  from  equations  (2). 
(3),  and  (6);  see  also  Medina  [1990]).  In  21-m  depth,  Middleton 
et  al.  [1987]  noted  a  similar  dependence  of  the  infragravity  energy 
on  the  sea  and  swell  ireciuencies.  Infragravity  energy  increased 
when  swell  (11-14  s)  energy  increased,  but  not  with  increases  in 
shorter-period  (SIO  s)  waves.  A  dqtendence  of  the  infragravity 
energy  levels  on  the  incident  sea  and  frcqu.'nwks  has  also 
been  observed  in  the  surf  zone  [Goda,  1975;  Holman  and  Sal- 
lenger,  1985]. 

In  deep  water  (kk  »  1)  for  collinear  free  waves. 


,G  (deep)  —  — -j-(A„*y  —  A*  )  —  2^^^*  XtttM  *-j  ~tt), )  (7) 

and  the  ratio  of  bound  infragravity  to  free  sea  and  swell  energy  is 
proportional  to  (akXaco^/g)  where  ak  is  the  first-order  free 
wave  steepness  and  Wy  is  the  bound  wave  frequency.  At  a  fixed 
®y,|,C|  increases  for  increasing  free  wave  frequencies  (or 
wavenumber)  in  deep  water,  and  sea  forces  more  infragravity 
energy  at  the  surface  than  is  forced  by  swell  of  the  same  height 
(compare  spectra  for  2(X)-m  depths  in  Figures  3a  and  3f>);  the  rela¬ 
tive  response  is  reversed  compared  to  shallow  water.  Owing  to 
the  relative  importance  of  steep  wind  waves  to  infragravity  wave 
generation  in  deep  water,  subsurface  pressure  sensors  (which  can¬ 
not  sense  short  wind  waves  because  of  hydrodynamic  attenuation) 
may  not  provide  data  suitable  for  bound  wave  predictions  at  the 
surface.  Also  in  contrast  to  shallow  water,  in  deep  water  |,C  | 
increases  with  increasing  bound  wave  frequency  so  the  spectrum 
of  bound  waves  in  deep  water  is  generally  not  white  but  most 
energetic  at  a  frequency  that  roughly  corresponds  to  the  width  of 
the  free  wave  spectrum. 

3.  Field  Measvrfmekis  And  Da  i  a  Pkocessino 

Field  measurements  were  made  at  Barbers  Point  Harbor, 
Hawaii,  on  the  southwest  side  of  Oahu;  at  Imperial  Beach,  Cali¬ 
fornia;  and  offshore  of  Point  Conception.  California.  The  Barbers 
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Fig.  2.  Bound  wave  sea  surface  'Jevatton  spectra  (dashed  lines)  and  bottom  (z  ~-h)  pressure  (dotted  lines)  in  depths  of  (a)  1 0  m, 
(h)  40  m,  (c)  183  m  and  (d)  5000  m  forced  by  the  same  first-order  free  wave  (solid  line,  /  >0,04  H/.)  energy  spectrum  (Figure  2a) 
and  directional  distribution  (all  frequency  bands  have  the  same  directional  spread  determined  by  d).  The  full  angular 

widths  (degrees)  at  50%  of  the  maximum  free  wave  spectral  energy  density  ate  indicated,  and  the  tout  sea  and  swell  energy  is  lO’ 
cm^.  The  '' '  '  bound  wave  sea  surface  elevation  and  bottom  pressure  energy  for  f  <0.04  Hz  for  the  various  angular  widths  and 
depths  are  given  in  Table  1. 


Point  data  were  collected  between  lune  1988  and  Match  1990  as 
part  of  a  study  of  oscillations  within  the  0.37  Icm^,  n.6-m-deep 
harbor  basin  (Figure  4).  Access  to  th^  basin  is  through  a  f.3-kin- 
long  entrance  channel,  dredged  to  a  depth  of  12.8  m.  The  coral 
rock  bearch  slopes  rather  steeply  (-0.04  slope)  from  the  shoreline 


to  3-ni  depth,  then  slopes  gradually  (~0.01)  to  the  9-jn  contour 
approximately  0.8  km  offshore,  followed  by  a  steeps-  (-0.08) 
slope  between  the  9-m  and  18-m  contours,  beyond  which  the 
depth  increases  to  180  m  within  another  0.8  km.  We  primarily 
use  data  from  four  near-bottom  pressure  sensors  arranged  in  a  6  m 
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TABLE  1.  Touil  Bound  Wave  Sea  Surface  Elevation  and  Bottom  Presture  Energy  for  /  <0.04  Hz  for  ihe  Angular  Widihs  and  Depihi  in  l  igure  2 


Depth,  m 

Angular  Width 

160° 

_ _ 

icr 

_ 40° _ 

Sea  Surface 
Elevation 
Energy,  cm^ 

Bottom 
Pressure 
Energy,  cm^ 

Sea  Surface 
Elevation 
Energy,  cm^ 

Bottom 
Preisure 
Energy,  cm^ 

Sea  Surface 
Elevation 
Energy,  cm^ 

Bottom 
Pressure 
Energy,  cm^ 

Sea  Surface 
Elevation 
Energy,  cm^ 

Bouom 
Pressure 
Energy,  cm^ 

10 

41.0 

40.6 

10.2 

10.0 

3.1 

3.0 

1.1 

Ll 

40 

0.16 

0.15 

0.075 

0.066 

0.028 

0,024 

0.031 

9(hlfr’ 

183 

0.0031 

3.7x10-^ 

0.0024 

1.8x10^ 

0.0025 

6,2x10’ 

0,016 

2.2*10’ 

5000 

0.0027 

<10-’ 

0.0020 

<10-’ 

0.0024 

0.016 

<10'' 

X  6  m  square  slope  array  located  outside  the  harbor  just  north  of 
the  entrance  channel  in  8.7  m  mean  depth.  Data  were  also 
acquired  from  bottom-mounted  pressure  sensors  in  7.0-  and  4.5-m 
depth  (Figure  4).  The  tidal  range  at  Barbers  Point  was  approxi¬ 
mately  0.9  m. 

Imperial  Beach,  located  just  north  of  the  Califomia-Mcxico 
border,  is  a  west-facing  sandy  beach  with  a  relatively  straight 
shoreline  and  plane  beach  slope  (-0.025).  The  Imperial  Beach 
data  set  is  from  a  6  m  x  6  m  slope  array  located  approximately 
0.4  km  offshore  in  about  11.5-m  mean  depth,  collected  from 
December  1989  to  August  1991.  The  tidal  range  was  approxi¬ 
mately  2.4  m.  At  Barbers  Point  and  Imperial  Beach  2.3 -hour-long 
records  were  obtained  at  a  1  Hz  sample  rate,  four  times  a  day. 
There  were  557  and  2301  records  at  Barbers  Point  and  Imperial 
Beach,  respectively. 

Offshore  data  were  obtained  with  a  Baylor  wave  staff  mounted 
on  the  Chevron  oil  platform  Hermosa  in  183-m  depth  approxi¬ 
mately  10  km  from  shore  between  Point  Arguello  and  Point  Con¬ 
ception,  California.  Ten  time  series  each  of  34-min  duration  and 
sampled  at  2  Hz  were  analyzed. 


The  energy  density  in  the  infragraviiy  band  can  be  10^  times 
lower  than  in  the  sea  and  swell  band  (c.g..  Figure  6),  and  'he  qual¬ 
ity  of  these  relatively  low  signals  could  be  degraded  by  nonlinear¬ 
ity  in  either  the  basic  sensor  (either  pressure  or  wave  staff)  or 
"flow  noise"  due  to  disturbance  of  the  flow  by  a  pressure  sensor. 
However,  infragravity  spectra  measured  by  pressure  sensors 
within  a  slope  array  were  similar,  as  would  probably  not  be  the 
case  if  sensor  nonlinearity  or  flow  noise  were  dominant 

Fourier  coefficients  of  pressure  and  sea  surface  elevation  were 
calculated  for  overlapped  1024-s  records  which  were  first  quadrat- 
ically  detrended  to  suppress  tidal  and  other  low-frequency 
motions  with  periods  greater  than  the  record  length,  and  tapered 
with  a  triangular  window  to  reduce  spectral  leakage.  The  Fourier 
coefficients  of  pressure  at  sea  and  swell  frequencies  were  con¬ 
verted  to  sea  surface  elevation  with  linear  theory.  Smoothed 
power  spectra  for  each  2.3-hour  record  were  obtained  by  averag¬ 
ing  over  1024-s  segments  and  merging  over  five  frequency  bands, 
resulting  in  a  spectral  bandwidth  of  4.88x10'^  Hz  with 
100-200  degrees  of  freedom  (DOF)  at  Barbers  Point  and  Imperial 
Beach  and  -30  IX)F  at  Point  Conception.  Frequency-directional 


(a)  (b) 


Fig.  3.  Bound  wave  sea  niiface  elevation  spectra  (dashed  lines)  and  botum  pressure  spectra  (dotted  lines)  forced  by  unidirectioiMl 
(a)  swell  (-0.07  Hz  peak)  and  (h)  sea  (-0.2  Hz  peak)  of  equal  toul  energy  (10^  cm^  solid  line)  in  depths  A  3.  10,  100,  and  200  m. 
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Fig.  4.  Ifiitnimem  posiUons  and  depths  at  Barbers  Point.  Solid  dots  are  single  point  pressure  sensors,  and  the  solid  dot  with  cro.ss 
is  a  6  m  X  6  m  array  of  four  pressure  sensors. 


spectra  were  crudely  estimated  at  Imperial  Beach  and  Baitiers 
Point.  The  mean  direction  Bo(/)  and  spreading  parameter  s(f) 
in  the  widely  used  cos^^'>l9(f)-0dlf)]  parametric  directional 
distribution  were  obtained  from  the  slope  array  cross  spectra  (see 
appendix). 

4.  Bound  Wave  PREDtenoNs  AND  Observations 

The  sea  surface  elevation  and  pressure  variance  below  0.04  Hz 
were  defined  as  inftagravity  wave  energy  and  , 

respectively),  while  the  sea  surface  elevation  variance  between 
0.04  and  0.25  Hz  was  assumed  to  be  free  sea  and  swell  wave 
energy  (Ess  )•  The  fixed  cutoff  frequency  of  0.04  Hz  between 
infragravity  and  sea  and  swell  waves  was  selected  so  that  very  lit¬ 
tle  energy  of  long  swell  was  included  in  the  infragravity  band.  A 
basic  question  that  is  addressed  is  whether  bound  waves  dominate 
the  total  energy  in  the  infragravity  band.  If  the  cutoff  frequency  is 
high  enough  so  that  significant  amounts  of  swell  energy  fall  in  the 
defined  infragravity  band,  then  the  bound  wave  prediction  will  be 
much  lower  than  observed,  and  the  result  will  necessarily  be  that 
bound  waves  are  insignificant  The  0.04-Hz  cutoff  frequency  was 
based  on  comparisons  of  many  predicted  bound  wave  and 
observed  spectra.  As  illustrated  in  Figures  6  and  12  (to  be  dis¬ 
cussed  later)  a  spectral  minimum  was  nearly  always  observed  at  a 
frequency  near  (or  above)  0.04  Hz.  At  frequencies  higher  than 
this  minimum,  the  observed  power  level  rises  dramatically  and  the 
predkied  bound  wave  energy  levels  are  much  lower  than 
observed  owing  to  the  importance  of  swell.  Obviously,  the  value 
of  a  fixed  cutoff  would  vary  according  to  the  wave  climate  at  a 
particular  site;  0.04  Hz  is  appropriate  for  the  present  observations 
in  the  Pacific  Ocesn. 


The  frequency  band  studied  here  (/  <0.04  Hz)  docs  not  neces¬ 
sarily  contain  most  of  the  energy  forced  by  nonlinear  difference 
interactions.  For  example,  in  100-  and  200-m  depths  in  Figure  3 
there  is  much  more  forced  energy  at  the  sea  surface  in  the  band 
0.04  to  0.20  Hz  than  below  0.04  Hz,  although  the  0.04  to  0.20-Hz 
forced  energy  is  much  lower  than  the  free  swell  energy  in  this 
band.  Thus  we  attempt  not  to  identify  all  the  energy  forced  by 
difference  frequetjcy  interactions,  but  rather  to  assess  whether 
forced  waves  dominate  a  low-frequency  band  where  free  waves 
might  be  negligible.  The  relatively  small  amount  of  bound  wave 
energy  above  0.04  Hz  (Figures  6  and  12)  included  as  part  of  the 
free  wave  (forcing)  spectrum  at  swell  frequencies  docs  not  cause 
significant  errors  in  bound  wave  predictions  below  0.04  Hz. 
There  are  cases  (not  shown)  in  which  the  spectrum  above  0.04  Hz 
consists  of  a  swell  peak  at  /,  and  an  additional  peak  at  2/p . 
Bispectral  analysis  [e.g.,  Elgar  and  Gaza,  1985]  shows  the  energy 
at  2/p  to  be  phase^oupled  bound  harmonics  driven  by  sum 
interactions,  and  not  free  wave  eitergy  as  assumed  in  the  calcula¬ 
tion  of  infragraviiy  bound  wave  energy.  This  is  of  little  conse¬ 
quence  here,  however,  because  in  shallow  water  the  swell  peak 
provides  virtually  all  the  bound  infragravity  wave  forcing  (Figures 
1  and  3)  and  the  predictions  are  insensitive  to  harmonic  energy. 

Similar  to  previous  results,  the  observed  infragravity  wave 
heights  (//'*•)  and  sea  and  swell  significant  heights  (Hss)  at  both 
Barbers  Point  (8.7 -m  depth)  and  Imperial  Beach  (10.2-12,7  m) 
are  significantly  correlated  (Table  2).  (Significant  heights  are 
defined  as  4<T,  where  <T*  is  the  band-passed  variance  in  each  fre¬ 
quency  range.)  H"**  and  Hss  are  also  highly  correlated  in  shal¬ 
lower  depths  (4.5  m  and  7.0  m)  at  Barbers  Point  and  deeper 
(183  m)  water  at  Point  Conception.  The  slope  of  the  best  fit  linear 
regression  line  through  the  data  points  (//<*•,// 55)  -O(0.i)  in 
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TABLE  2.  Simple  Relationships  Between  //”**  and  liss 


Siie 

DcfHh,  m 

Slope 

Con 

T.P./S.B 

swash 

1,00 

0.78 

T.P./S.B. 

0.0-  I.O 

0.37 

0.54 

T.P/S.B. 

O 

o 

032 

0.86 

B.P 

4.5 

0.20 

0.96 

B.P 

7.0 

0.17 

0.90 

B.P 

8.2 -9.1 

0.14 

0.92 

i.B. 

10.2-11.0 

0.13 

0.87 

LB. 

11.1  -  11.8 

0.13 

0.88 

LB. 

11.9-12.7 

0.13 

0.85 

LB. 

10.3-12.7 

013 

0.86 

North  Sea 

40.0 

0.04 

i>.C. 

183 

0.03 

0.97 

Slope  is  the  slope  of  the  linear  regression  line  between  //°**  and  Ilss,  and 
CoTT  is  the  correlation.  All  correlations  shown  are  significant.  Sites  in  the 
present  study  arc  Barbers  Point  (B.P.),  Imperial  Beach  G-B.),  and  Point 
Conception  (P.C.).  Similarly  processed  data  from  forrey  Pines  (T.P.)  and 
Santa  Barbara  (S.B.)  are  included  for  comparison;  Uss  in  these  cases  was 
measured  in  7- to  10-m  depth  (Giaa  and  Thornton,  1985J.  Sand  (I982i| 
had  two  records,  each  1  hour  long,  from  the  North  Sea. 

depths  of  -8-13  m  is  similar  to  past  observations  in  these  depths 
[e.g..  Munk,  1949;  Tucker,  1950],  as  are  the  shallower  observa¬ 
tions  at  Barbers  Point  (Table  2,  and  data  from  previous  studies  not 
shown).  Regression  line  slopes  are  considerably  smaller  in  both 
40-m  [Sand,  1982h]  and  183-m  depth.  Note  that  bound  wave 
theory  does  not  predict  a  linear  relatioriship  between  //"**  and 
ttss,  for  a  free  wave  field  with  the  shape  of  the  frequency- 
directional  spectrum  and  depth  held  constant,  and  only  the  total 
power  level  varied,  H'*‘  -  Hss  *• 

Wave  spectra  in  -10-m  depth  at  Imperial  Beach  and  Barbers 
Point  were  separated  into  three  groups  delineated  by  the  fre¬ 
quency  of  the  maximum  spectral  level.  As  shown  in  Figure  5,  for 
a  given  toul  sea  and  swell  energy  {Ess ).  the  observed  infragravity 
energy  (fE"**)  is  clearly  larger  for  swell  than  for  higher  frequency 
[(7(0.10  Hz)]  wind  waves  (particularly  for  energetic  events),  a 
trend  qualitatively  consistent  with  both  bound  wave  theory  (equa¬ 
tion  (6)  and  Figure  3)  and  other  observations  [Middleton  et  ai. 
1987;  Nelson  et  ai,  1988].  With  similar  Ess,  in  183-m 
depth  is  near  10'^  times  less  than  in  10-  to  13-m  depth  (Figure  Sb). 

The  observed  sea  surface  elevation  spectra  above  0.04  Hz 
(assumed  to  be  free  waves)  were  used  to  predict  (equation  (3))  the 
infragravity  bound  wave  spectra.  Figure  6  (and  the  top  panel  of 
Figure  12,  discussed  later)  shows  examples  of  observed  and 
predicted  pressure  spectra  with  the  sea  and  swell  directional  spec¬ 
tra  modeled  both  as  cas^^^\9{f )  -  floC/" )]  (appendix)  and  as  uni¬ 
directional. 

Owing  to  the  maximum  of  the  bound  wave  coupling  coefficient 
C  for  collinear  waves  (equation  (3),  Figure  2a),  pE*^  is 
between  2  and  10  times  larger  with  unidirectional  free  waves  (sea 
and  swell)  than  with  directionally  spread  free  waves.  Bound 
waves  were  also  predicted  assuming  the  sea  and  swell  directions 
were  normally  distributed  about  flo  with  rms  directional  width 
given  by  equation  (A2).  Differences  between  the  predictions 
using  the  cos^  and  Gaussian  directional  distribution  were  small. 
However,  the  sensitivity  of  the  bound  wave  model  to  the  direc¬ 
tional  distribution  is  strong  (Figure  2<i),  so  the  fundamentally  low 
resolution  of  the  slope  array  limits  the  accuracy  of  bound 
infragravity  wave  predictions.  Thus  the  unidirectional  wave 
pE^  is  an  upper  limit  on  bound  wave  model  predictions,  and 


^£t>nd  j  cos^  distribution  is  a  rough  approximation  of  the 
effects  of  directional  spreading. 

The  toul  observed  (/  <0.04  Hz)  infragravity  energy  ipE'*^) 
for  the  entire  Barbers  Point  and  Imperial  Beach  daU  sets  arc  com¬ 
pared  with  bound  wave  predictions  {pE^)  with  both  unidirec¬ 
tional  and  directionally  spread  free  waves  in  Figure  7.  With  direc¬ 
tionally  spread  waves  (Figures  la  and  lb),  bound  wave  theory 
always  predicts  less  energy  than  observed  {pE'^  >  pE^).  This  is 
not  always  so  with  unidirectional  free  waves  (Figures  7c  and  Id), 
emphasizing  the  importance  of  including  directional  spreading  in 
bound  wave  calculations. 

Both  directional  and  unidirectional  models  severely  undcr- 
prcdict  {pE"^  »  pE'*^)  when  pE‘*'  (and  sea  and  swell 

wave  energy;  Figure  5)  is  small.  The  underprediction  when  uni¬ 
directional  sea  and  swell  waves  are  assumed  (which  results  in  an 
overestimate  of  pE^)  shows  that  other  types  of  infragravity 
motions  besides  bound  waves  must  sometimes  be  important  in 
these  data  sets.  The  following  discussion  considers  pE*^ 
predicted  with  directionally  spread  waves  only,  keeping  in  mirxl 
the  low  resolution  of  the  directional  estimates  and  the  sensitivity 
of  bound  wave  theory  to  directional  spreading. 

Much  of  the  scatter  in  ,£°*”  for  fixed  pE'^  at  Imperial  Beach 
(Figure  la)  appears  to  be  correlated  with  depth  changes.  The  tidal 
range  is  2.4  m  at  Imperial  Beach  compared  with  0.9  m  at  Barbers 
Point,  and  there  is  indeed  less  scatter  at  Barbers  Point  (compare 
Figures  la  and  lb).  The  underprediction  of  at  Imperial 

Beach  is  most  pronounced  at  high  tide  (depth  of  12-13  m)  with 
small  swell  and  p£°**,  and  least  at  low  tide  with  energetic  p£** 
(spectra  for  representative  extreme  cases  are  shown  in  Figure  6). 
The  data  were  binned  according  to  mean  water  depth  and  pE’**. 
The  ratios  (averaged  within  a  bin)  of  predicted  bound  wave 
energy  to  observed  infragravity  energy  (,£‘™‘/p£'**),  an  estimate 
of  the  fraction  of  the  total  infragravity  energy  which  is  bound  in 
each  bia  increase  nearly  monotonically  with  increasing  pE‘*‘ 
and  decreasing  depth  (Figure  8,  Table  3).  On  average,  less  than 
5%  of  the  energy  is  estimated  to  be  bound  when  the  depth  h  is 
greater  than  11.5  m  and  pE*’ <6  cm^,  but  more  than  20%  is 
bound  when  h  <11  m  and  >10^.  The  similar  values  of 

pE'^IpE*’  in  the  deepest  depth  bin  at  Barbers  Point 
(850-913  cm)  and  the  shallowest  Imperial  Beach  depth  bin 
(100&-1050  cm)  are  noteworthy  given  the  differences  between  the 
sites  (i.e.,  the  mainland  coast  at  Imperial  Beach  compared  with  the 
island  site  near  a  harbor  at  Barbers  Point).  Because  pE’^  and 
Ess  ere  strongly  correlated  (e.g.,p£”’”  -10  cm^  when  Ess  -ICXX) 
cm^;  Figure  5),  Figure  8  implies  that  the  fraction  of  the  infragrav¬ 
ity  energy  which  is  bound  (pE’^/pE'*')  also  increases  with 
increasing  Ess  -  This  is  consistent  with  the  results  of  Elgar  el  ai 
[1992]  in  8-  to  13-m  water  depth  at  Duck,  Nwth  Carolina. 

Although  the  ratios  averaged  over  bins  show  smooth  trends 
with  depth  and  pE**  (Figure  8),  there  is  considerable  scatter  of 
the  ratios  within  a  bin  (Figure  9).  Additional  variables  besides 
depth  and  p£“*»  apparently  influence  the  fraction  of  the  energy 
which  is  bound.  For  example,  certain  directional  characteristics 
of  the  incident  wave  field  may  preferentially  excite  edge  waves 
[Gallagher,  1971]  and  thus  reduce  the  relative  importance  of 
bound  waves.  Note  that  the  tendency  for  higher  pE"**  to  occur 
with  low-frequency  incident  swell  than  with  higher-frequency 
seas  of  the  same  energy  (Figure  5)  is  not  apparent  in  the  ratio  data 
(Figure  9).  Thus  the  fraction  of  the  energy  which  is  bmind  does 
not  show  a  clear  trend  with  the  peak  frequency  of  the  incident 
waves. 

At  the  offshore  Point  Conception  site  in  183-m  depth,  there 
was  only  a  single  wave  gage,  so  only  unidirectitmal  model  predic- 
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Fig.  5.  Observed  infragnvity  (0.005  </  <  0.04  Hz)  pressuie  energy  (pE**)  venui  tei  end  swell  (0.04  </  <  0.25  Hz)  surface 
elevation  energy  (Mss)  st  (a)  Barbers  P^t,  8.2  to  9.t-m  depth;  (b)  Imperial  Beach,  10.2  to  12.7-m  depth.  Each  plotted  point  is 
from  a  2.3-haur  record.  The  dashed  line  it  pE^^mEss  where  ns  is  the  slope  of  the  best  fit  linear  regression  line  between 
pE"**  and  Ess  i  "»  •  0.016  (0.017)  at  Barbers  Point  Ompetial  Beach),  and  the  conelation  between  p  E*"  aird  Ess  ts  0.88  (0.83). 
Peak  ftequenciet  (fp)  of  each  record  are  /,  S0.067  Hz  (asterisks),  0.067  <  /,  SO.IO  Hz  (circles),  and  fp  >0.10  Hz  (triangles). 
The  10  Point  Conception  obseivatiaiu  (183-m  depth)  of  and  Ess  »•*  shown  by  crossed  squares  (Figure  5b). 


lions  were  made  (Figure  10).  The  average  (ova-  all  10  records) 
ratio  is  0.047,  and  accounting  for  the  directional 

spread  in  the  sea  and  swell  further  decreases  the  energy  ratio 
(unless  the  directional  distribution  is  unusually  broad  (Figure  2c, 
Table  1)).  The  incklent  sea  and  swell  were  energetic  at  Point 
Conception,  relative  to  typical  conditions  at  Imperial  Beach  and 


Barbers  Point  At  Point  Conception,  Ess  >1400  cm^  for  all  10 
records,  and  the  average  Ess  -  4400  cm^  which  is  comparable  to 
the  largest  Ess  u  the  shallow  sites  (Figure  S).  Thus  on  average 
only  a  small  amount  (<5%)  of  .£“•■  is  bound  in  183-m  depth, 
even  with  energetic  incident  waves  for  which  >25%  of  the  energy 
is  bound  in  shallow  water  (FigtHcs  S  and  8). 
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Fig.  6.  Observed  see  and  swell  surface  elevation  spectnim  (solid  line,  /  >  0.04  Hr.),  observed  bollcm  pressure  spcctnim  (solid 
line,  /  <0.04  Hz),  unidirectional  bound  wave  model  predicted  pressure  spectrum  (dashed  line),  and  directional  (cos^)  model 
predicted  pressure  spectrum  (dotted  line)  at  Imperial  Beach,  (a)  March  3.  1991,  0710,  h  =  10.9  m.  Ess  =  2239  cm^  =  71 
cm^  and  f,E^  =110  and  28  cm^  for  unidirectional  and  directional  mt  'el,  respectively,  (fr)  November  14,  1990,  0702;  A  =  12.5 
m.  Ess  =  120  cm^,  ,£*'=  1.6  cm^  and  ,,£'*“*=  0.09  and  0.009  cm^  for  unidirectional  and  directional  model,  respectively. 
The  pE°**  and  p£*"**  ate  based  on  /  <  0.04  Hz. 


1 


5.  Discussion 

Wave  Type  Mix 

The  predicted  bound  wave  energy  (£*^)  is  less  than  the 
observed  total  infragravity  energy  (£'’’*),  both  in  shallow  (~10-m 
depth.  Figures  7a  and  76)  and  deep  (183-m  depth.  Figure  10) 
water.  The  additional  observed  energy  is  contributed  by  free 
waves  (£**),  consisting  of  both  leaky  (£''’^)  and  edge  (£“*») 
waves. 

£<*•  =  £'«'  +  £«”,  (8a) 

where 

£f»>  =  £iJ>r  +  £«4a.  (86) 

By  comparing  £•*•  in  different  depths  (6i  and  62,6)  <62),  for 
similar  sea  and  swell  wave  energy  (Ess),  the  fraction  (Li,  1,1,2)  of 


ftrs  in  depth  6|  which  radiates  to  a  deeper  site  (depth  62)  can  be 
estimated.  The  assumption  is  made  that  Ess  and  h  determine 
the  total  infragravity  energy  (and  mis  of  wave  types),  so  observa¬ 
tions  at  the  three  sites  may  be  combined  to  form  an  incoherent 
array.  £”•*  and  the  average  fraction  of  £“*•  estimated  to  be 
bound  (B ), 

£b»d  =  g£,M  =  _  (9) 

arc  given  in  Table  4  for  the  subset  of  the  data  with  Ess  about 
1900  cm^  (Hss  ~175  cm).  The  small  number  of  observations  in 
deep  water  (Point  Conception,  Figure  56)  is  limited  to  cases  of 
energetic  sea  and  swell. 

It  follows  from  the  conservation  of  energy  (EC,  =  constant, 
where  C,  is  the  group  velocity)  that  leaky  infragravity  energy 
E'^  ~  6"“  for  normally  incident  waves  in  shallow  (relative  to  the 
infragravity  wavelength)  water.  On  plane  parallel  contours. 
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Fig.  7.  (a,  6)  DiicctiofMl  and  (c,  d)  unidireciicntl  bound  wive  model  predicted  inftagravity  presiure  energy  versus 

obterved  infragnvity  preisure  energy  (^£°^)  at  (a,  c)  Imperial  Beach  and  (h,  d)  Barbers  Point.  The  symbols  indicate  the  mean 
depth  (A  )of  each  record:  diamonds,  820  ^  S913  cm;  K)uares,  1000  ^  $1099  cm;  pluses,  1100  SA  $1199  cm;  mangles,  1200 
$A  $1300  cm.  The  solid  lines  are  =  ,£■*•.  To  improve  figure  darity,  only  one  fourth  of  the  data  poinu  at  Imperial  Beach 
are  plotted,  so  about  600  points  are  shown  for  each  site. 


moderate  angles  of  oUiqoity  (<4S°)  in  deep  water  cause  less  than 
30%  deviation  from  the  k  dependence.  Thus 

£»7-(Aa)  =  £M/it)[/>i/A2r.  (10) 

The  refractive  trapping  of  edge  wave  energy  between  depths  A, 
and  Ai  will  reduce  the  amount  of  £*"  reaching  Aj  by  the  factor 

£»«(Ai)  =  £«-(Ai)(A  i/Ail'*£»,,*. .  (11) 


Manipulation  of  equations  (8a),  (9),  and  (1 1)  yields 


1-«(A2) 

£'**(A2) 

hz 

1  ~B  (Ai) 

£**(Ai) 

(12) 


When  B  is  small,  as  in  the  present  case  (Table  4),  Lh^  ,a^  is  insen¬ 
sitive  to  large  relative  changes  in  B  and  is  essentially  a  measure 
of  the  deviation  of  £*•  from  an  A"'*  leaky  wave  dependence. 
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Fig.  S.  Coniourj  of  Ihe  ratio  of  predicted  bound  wave  energy  to  observed  infragravily  energy  as  a  function  of  depth 

and  pE**,  at  (top)  Barbers  Point  and  (bottom)  Imperial  Beach.  Data  were  binned  according  to  mean  water  depth  (eight  bins)  and 
observed  infragravity  energy  (five  logarithmically  spaced  bins)  and  ratios  averaged  within  a  bin.  The  average  ratios  and  number  of 
observations  in  each  bin  are  given  in  Table  3. 


Substitution  of  observed  values  (Table  4)  in  equation  (12)  yields 
I'tJ.ni—  1.02,  Laj.iw  =1-11.5, 1S3  =  0.07.  Virtually  no  free 
wave  energy  is  trapped  between  8.5-  and  11.5-m  depth,  compared 
with  greater  than  90%  between  either  shallow  sensor  and  183-m 
depth.  Thus  in  8-  to  13-m  depth,  for  Ess  -1900  cm^,  an  average 
of  about  25%  of  the  total  infragravity  energy  is  estimated  to  be 
bound  (Table  4);  -5%  is  leaky  (or  of  high  mode  number)  and 
reaches  183  m  depth,  and  70%  is  trapped  shoreward  of  183  m. 
Frequency- wavenumber  spectra  from  an  array  of  pressure  sensors 
in  8-m  depth  at  Duck,  North  Carolina,  also  suggest  die  importance 


of  edge  waves  well  outside  the  surf  zone  [OUman-Shay  et  ai, 
1989;  J.  Oltman-Shay,  personal  communication,  1990]. 

The  spatial  variations  of  and  £“**  for  idealized 

wave  conditions  are  illustrated  in  Figure  11.  E^  is  calculated 
using  equation  (3)  and  the  frequency  spectrum  shown  in  Figure  2 
with  Ess  =  1800 cm^  (close  to  Ess  in  Table  4)  and  full  angular 
width  of  40”  (selected  so  the  predicted  and  observed  £*^  in 
1 1.5-m  depth  were  approximately  6  cm^;  ££<*•  in  Table  4).  E^ 
for  a  typical  infragravity  frequency  (0.02  Hz)  and  normal 
incidence  i.s  also  shown  in  Figure  11.  Shoreward  of  the  turning 


TABLE  3.  Average  Ratios  pE^tpE'**  and  Number  of  Observations  in  Each  Bin  in  Figure  8 

___ 


Depth,  cm 

0.0-1 .0 

1.0-3.16 

3.16-10.0 

10.0-31.6 

31.6-100.0 

Ratio 

No. 

Ratio 

No. 

Ratio 

No. 

Ratio 

No. 

Ratio 

No. 

820-849 

0.09 

71 

0.15 

73 

0.22 

20 

0.25 

14 

0.26 

3 

850-913 

0.08 

117 

0.13 

187 

0.21 

42 

0.22 

25 

0.36 

5 

1000-1049 

0.07 

9 

0.13 

13 

0.23 

6 

0.46 

2 

* 

1050-1099 

0.05 

87 

0.10 

146 

0.18 

69 

0.31 

38 

0.35 

4 

1100-1149 

0.03 

112 

0.04 

307 

0.08 

213 

0.17 

91 

0.31 

22 

115(7-1199 

0.01 

40 

0.02 

352 

0»M 

305 

0.11 

84 

0.20 

20 

1200-1249 

0.01 

6 

O.Ol 

128 

0.02 

162 

007 

48 

0.16 

7 

1250-1300 

0.01 

15 

0.02 

9 

0.06 

5 

0.16 
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Fig.  9.  The  niio  of  pR<l)aed  bound  wm  to  obierved  tnfngraviiy  energy  (pE^pE'*’),  (a)  as  a  function  of  pE'*’  for  the  lim¬ 
ited  depth  range  1050  -  1 100  cm  and  (i)  as  a  function  of  depth  for  a  fixed  range  of  pE’*‘  (10-32  cm^).  At  in  Figure  5.  peak  fre¬ 
quencies  {fp)a(  each  record  are  fp  S0.067  IIz  (attehski),  0.067  <fp  S0.10  Hz  (circlet)  and  fp  >0  10  I!/,  (triangles).  The  dathed 
vertical  lines  indicate  bin  boundaries  for  tite  average  ratios  shown  in  Figure  8. 


point  (the  most  seaward  edge  wave  antinode  where  exponential 
decay  begins),  and  beyond  the  turning  point, 

E«da  _  g-’V'* ,  where  k,  is  the  longshore  wavenumber  and  x  is 
the  aoss-shore  distance  measured  from  the  shoreline.  Thus  leaky 
and  edge  waves  have  simitar  spatial  dependence  shoreward  of  me 
edge  wave  turning  point.  The  exponentially  decaying  edge  wave 


Fig.  10.  Unidirectional  bound  wave  ntodel  predicted  infragravity  aea  sur¬ 
face  elevation  energy  (pE*^  vetius  obterved  infragravity  tea  surface 
elevation  energy  (,£***)  in  l83-m  depth  at  Point  Conception.  The  lolid 
line  is 


tails  in  Figure  1 1  show  that  of  the  edge  waves  modes  high  enough 
(n  i3)  to  reach  8.5-m  depth,  only  the  amplitude  of  a  mode  3  edge 
wave  varies  significantly  over  the  depth  range  8,5-1 1 ,5  m  (for  a 
frequency  0.02  Hz  on  a  0,02  slope).  Mode  numbers  fi7  are 
trapped  shoreward  of  183  m.  The  location  and  depth  of  a  turning 
point  depends  on  the  edge  wave  mode  number,  frequency,  and  the 
beach  slope  but  it  is  generally  the  case  that  many  more  modes 
have  turning  points  between  1 1.5  and  183  m  than  between  8.5  and 
Il.5mdepihs.  £*•  in  1 1.5-m  depth  (Figure  1 1 )  is  set  equal  to  the 
observed  value  (E”**  (l-B)  =  22  cm^  from  Table  4).  The  curves 
for  and  are  thus  both  constructed  to  match  the  observa¬ 
tions  in  11.5-m  depth.  The  calculated  Lts,M.s=1.0  (equation  (12) 
attd  Table  4)  suggests  little  trapping  of  infragravity  energy 
between  depths  of  8.5  and  11.5  m.  On  the  other  hand,  £*•  at 
Point  Conception  is  much  larger  than  would  occur  with  bound 
waves  only,  but  much  smaller  than  if  alt  the  free  wave  energy  in 
shallow  water  radiated  offshore.  The  calculated 
£ij,iai=£iu.m  =  0'07  essentially  corresponds  to  the  roduciion 
of  the  Point  Conception  observation  below  the  free  wave  curve  in 
Figure  11.  Although  leaky  (and  very  high -mode  (n  >7)  edge) 
wave  energy  is  only  a  small  fraction  of  £■*•  in  sh.  ’low  water,  it  is 
sufficient  to  completely  dominate  the  infragravity  wave  energy  in 
183-m  depth. 

Sand’s  (1982l>]  sea  surface  elevation  observations  in  40-m 
depth,  yielding  £^/£®^-1.0,  are  much  different  from  the 


TABLE  4.  Avenge  Obierved  Infragravity  Energy  (£**•)  and  Average  Bound  Fraction  B  =  (£'^/£‘**) 
for  Similar  Average  Sea  Swell  Energy  (Ess  ) _ _ 


Site 

Mon  Depth 
m 

No.  of 
Reooidi 

cm* 

<£*•) 

cm* 

B.P, 

8.3 

91 

1950 

35 

0.28 

IB, 

11.5 

13 

1830 

28 

0.21 

PC. 

U3 

4 

1840 

0.43 

0.(77 

D«U  from  each  eitr  tie  averaged  together.  Silee  are  identified  in  the  focttiote  to  Table  2. 
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Fig.  11.  Infragrivity  w»ve  forms.  The  solid  line  is  the  bound  w»ve  energy  forced  by  i  sea  and  swell  spectrum  =  ISOOcm^ 
and  fp  =  0.067  Ha)  with  a  constant  40°  angular  width.  Leaky  wave  and  edge  wave  energy  envelopes  at  0.02  lU  on  a  boaom 
slope  of  0.02  are  shown  by  doaed  lines.  The  leakyxdge  wave  curve  it  constructed  such  that  the  ratio  of  bound  wave  energy  to  the 
total  Ocakytedge  +  bound)  infragravity  energy  is  -0.2  in  ll.S  m.  The  leaky.edge  wave  energy  at  other  depths  is  obuin^  from 
linear  conservation  of  energy.  L^y  and  edge  wave  curves  are  similar  shoreward  of  the  exponentially  decaying  edge  wave  uils, 
shown  for  modes  2.  ‘iJ5,  and  7.  The  asterisks  indicate  the  average  observed  infragravity  energy  (£■*•)  at  each  site,  for  sea  and 
swell  spectra  with  Ess  -1900  cm*  (Table  4).  The  depths  discussed  in  the  text  arc  shown  by  vertical  da  '  cd  lines. 


present  ratios  showing  a  maximum  value  of  0,15  in  12.5-m  depth 
and  a  decrease  with  increasing  depth  (Figure  8).  This  discrepancy 
may  be  due  to  a  number  of  factors.  First,  Sand's  predicted 
was  calculated  assuming  unidirectional  sea  and  swell  waves, 
which  overestimates  arxl  Although  the  very 

sti®ng  effects  of  directional  spreading  are  reduced  in  deep  water 
compared  with  shallow  water,  an  angular  .spread  of  40“  reduces 
the  predicted  E*^  by  about  a  factor  of  6  (Table  1).  This  is 
significant  and  may  explain  the  large  differences  (at  least  a  factor 
of  5)  with  the  present  observations.  However,  differences  in  the 
sites  and  Ess  are  probably  also  important.  Sand’s  observations 
were  from  the  middle  of  the  North  Sea.  about  200  km  from  the 
nearest  coastline,  presumably  reached  only  by  leaky  waves.  Thus 
edge  waves  present  at  the  current  sites  (on  the  continental  shelf 
and  in  Hawaii)  may  not  have  occurred  in  Sand's  observations,  and 
bound  waves  were  thus  relatively  more  important.  Additionally, 
E*^/E°*“  increases  with  increasing  Ess  (at  least  in  shallow 
water,  Figure  8)  and  Ess  -26,000  cm^  (Hss  -650  cm!)  in  the 
North  Sea  data,  coiisiderably  larger  than  the  maximum  Ess  -50(X) 
cm^  in  the  present  data  sets  (Figure  5).  Negligible  £“■«  and  large 
Ess  may  have  increased  the  relative  importance  of  bound  waves 
in  the  North  Sea. 

At  5000-m  depth  in  the  North  Atlantic  -500  km  from  the 
nearest  coastline,  Webb  el  aL  [1991]  found  that  E'*"  on  the 
seafloor  was  correlated  not  with  the  local  Ess  but  with  Ess  aver¬ 
aged  over  all  coastlines  within  the  'Tine  of  sight"  of  the  experi¬ 
ment  site.  The  conelations  occur  because  locally  forced  bound 
waves  are  too  short  to  penetrate  the  water  column  and  ,  E**  is 
due  predominantly  to  free  inffagravity  energy  £•",  radiated  from 
sometimes  distant  coastal  areas  with  energetic  sea  and  swell 


[Webb  el  al.,  1991).  In  contrast,  infragravity  waves  measured  on 
the  seafloor  in  40()0-m  depths  200  km  offshore  of  Point  Arena, 
California  [Sullon  el  al.,  1965),  and  at  the  sea  surface  in  183-m 
depth  10  km  west  of  Point  Concepcion  were  correlated  with  the 
local  Ess  even  though  bound  waves  are  not  important  in  either 
data  set.  E'*“  and  local  Ess  are  probably  correlated  because  the 
local  Ess  is  highly  correlated  with  Ess  at  nearby  coastal  sites 
radiating  free  infragravity  energy,  and  the  decrease  in  £^  due  to 
radial  spreading  is  minimal  because  of  the  proximity  of  these  sites 
to  Che  coast.  In  SOOO-m  depth,  lOOO  km  from  the  coast  and  due 
west  of  the  183-m  depth  site  of  the  present  study,  directional  spec¬ 
tra  [Webb  el  al.,  1991]  do  in  fact  show  significant  infragravity 
energy  propagating  from  the  California  coast  near  Point  Concep¬ 
tion,  where  surface  wave  energy  is  consistently  high.  The  present 
Point  Conception  data  set  contains  no  cases  of  small  Ess ,  so  radi¬ 
ation  from  the  nearby  coast  is  presumably  always  significant  and 
more  important  than  radia'ion  from  distant  coasts.  In  both  the 
North  Atlantic  (with  Ess  a  line  of  sight  average,  Webb  et  al., 
1991)  and  offshore  of  Point  Conception  in  the  Pacific,  a  very 
small  fraction  of  Ess  is  radiated  from  coastal  regions  as  leaky 
(and/or  very  high  mode  edge)  infragravity  energy. 


Third  Moments 

Theoretically,  infragravity  bound  waves  generated  by  nonlinear 
difference  interacu'ons  are  phase  coupled  to  pairs  of  free  waves 
(sea  and  swell)  and  180“  out  of  phase  with  the  wave  group 
envelope  (for  narrow  free  wave  directional  distributions;  negative 
C  in  Figure  1).  This  coupling  produces  negative  skewness. 
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TABLE  5.  Average  SkewneK  Values  for  34-min-I.ong  Recordi 


Sile 

Depth,  m 

No.  <rf 

Records 

Skewness 

1,5/, 

Nyquist 

B.P. 

8.2 -9.1 

5S7  06) 

43.08  (4).15) 

0.12  (028) 

I.B. 

10.2-11.0 

479  (103) 

-0.09  (-0.15) 

0.02  (0,04) 

LB. 

11.1  -  11.8 

1313(233) 

-0.07  (-0,14) 

0.02  (0.03) 

I.B. 

11.9-  12.7 

509  (82) 

-0.06  (-0.14) 

0,01  (003) 

RC. 

183 

10(10) 

0,003  (0.003) 

0.08  (0.08) 

The  upper  limits  of  frequencies  included  in  the  calcuUiions  are  indicated  in  the  heading  of  the  skewness 
columns.  The  Nyquist  frequency  was  1 .0  Hz  at  Point  Canception  (P.C)  and  0.5  Hz  al  Baihcn  Point  (B.P.)  and 
Imperial  Beach  (I.B.).  Energy  at  frequencies  greater  than  the  upper  limit  was  set  to  zero.  Values  in  parentheses 
are  for  the  larger  waves  (Ess  S62S  cm). 


Coupling  between  sea  and  swell  and  higher-frequency  waves 
forced  by  sum  interactions  produces  positive  skewness  [Hassel- 
mannetai,  1963). 

Data  from  all  three  field  sites  were  divided  into  34-min  seg¬ 
ments.  To  reduce  the  effects  of  high-frequency  bound  harmonics 
on  skewness,  each  record  was  low-pass  filtered  by  setting  equal  to 
zero  the  Fourier  coefficients  of  sea  surface  elevation  at  firequen- 
cies  greater  than  1.S  times  the  power  spectral  peak.  After  an 
inverse  Fourier  transform,  the  skewness  was  calculated  in  the  time 
domain.  Skewness  was  also  calculated  without  low-pass  filtering, 
thus  including  both  high-  and  low-frequency  bound  waves.  As 
shown  in  Table  S,  the  skewness  of  the  low-passed  data  is  negative 
at  the  shallow  depths  (8-13  m).  consistent  with  the  theoretically 
expected  coupling  between  sea  and  swell  and  infragraviiy  bound 
waves,  and  close  to  zero  in  deep  water  (183'm  depth)  as  expected 
for  uncoupled  infiragravity  waves.  The  magnitude  of  the  skewness 
in  shallow  water  increases  when  the  low-energy  records  (Ess  < 
625  cm^)  are  excluded  (Table  5).  Although  the  near-zero  band- 
passed  skewness  in  183-m  depth  is  consistent  with  negligible 
it  must  be  interpreted  cautiously.  In  this  depth, 
<!>■  (A  /g  ~  2,  directional  spreads  may  be  large,  and  bound 
infiragravity  waves  may  contribute  either  positive  or  negative 
skewness  (Figure  1).  Zero  skewness  could  thus  also  result  horn 
canceling  contributions. 

In  shallow  water  the  band-passed  skewness  (1 .5  column  in 
Table  5)  is  negative  while  the  total  (Nyquist)  skewness  is  positive. 
The  positive  skewness  due  to  bound  high-frequency  harmonics  is 
significantly  offset  by  negative  skewness  from  bound  infragravity 
waves,  similar  to  results  from  the  surf  zone  [Elgar  and  Guza, 
1985).  The  cancellation  is  nearly  complete  at  Imperial  Beach, 
where  positive  and  negative  contributions  of  about  |0.I  |  yield  a 
total  skewness  of  roughly  0.03. 

Bispectral  analysis  [Hasselmann  et  al.,  1963]  indicates  the  con¬ 
tribution  to  third  moments  from  indtvidual  triads  of  infragraviiy 
and  sea  and  swell  fiequencies.  The  normalized  magnitude  of  the 
bispectrum,  die  bicoherence  (h),  indicates  the  relative  degree  of 
phase  coupling  between  fiequencies.  The  biphase  gives  the  phase 
relationship  between  the  coupled  waves.  The  same  34-min 
records  for  which  skewness  was  calculated  were  analyzed  using 
bispectral  techniques.  An  example  from  Imperial  Beach  (Figure 
12)  shows  that  infiragravity  energy  is  significantly  coupled  to 
higher-finquency  sea  and  swell  within  the  power  spectral  pedt 
(0.05  S/  S  0.09  Hz).  The  biphases  (not  shown)  are  between 
158*  and  176*.  Rgtire  12  is  typical  of  cases  with  energetic  swell 
and  infragravity  waves  in  relatively  shallow  water  (4-  to  13-m 
depth). 

In  relatively  deep  water  (183  m;  ofishoie  of  Point  Canception) 
bispectra  do  not  ahow  statistically  significant  coiqiling  between 
sea  md  swell  and  infiragravity  waves,  consistent  with  our  result 
that  bound  waves  make  a  very  small  contribution  to  the  observed 
infiragravity  wave  spectrum. 


6.  Summary 

The  seiKitivity  of  the  predicted  [Hasselmann,  1962)  infragrav- 
ily  bound  wave  energy  (E“)  to  variations  in  d^lh  and  sea  and 
swell  conditions  was  examined  and  E^  was  compared  with  the 
observed  infragravity  energy  (£“**)  in  both  8-  to  13-m  and  l83-m 
water  depths. 

E^  varies  considerably  wiih  changes  in  the  depth,  the  free 
wave  (sea  and  swell)  frequencies,  and  free  wave  angular  spread. 
With  the  same  sea  and  swell  spectrum,  in  shallow  (0(10 
m))  and  deep  (O  (200  m))  water  can  differ  by  10*  (Figures  2  and 


Fig.  12.  (top)  Energy  spectra  and  (boaoro)  contoun  of  bicoherence  (24 
DOF)  in  10.9-m  depth  at  Imperil  BeaiA  (Januaiy  14,  1990,  1948). 
Observed  sea  and  swell  surface  elevation  spectnim  (solid  line,  /  i  0.04 
Hz),  observed  infngravity  pressure  spectrum  (stdid  In*.  /  <0,04  Hz),  uni¬ 
directional  bound  wave  mc^  preasure  spectium  (dashed  line),  and  cos^ 
model  pressure  spectrum  (dotted  line)  are  shown  (Ess  “  1^1^ 

27  cm\  31  and  13  cm*  for  unidireaional  and  diteakmal 

model,  respectively).  The  minimian  bicoheteiKe-contoor  level  (rfotied  is 
significant  at  the  95%  level  lHaiibr,x.:t,  1965],  b  e  0.S0,  with  additional 
contours  every  0.1.  The  convention  is  that  the  interacting  frequencies  are 
/i>  /a*  / 1  +/ 1,  with  fi  2/2:  thus  the  lowest  frequency  is  fi  and  the 
bicoheience  eiKlosed  in  the  dashed  box  indicate  coupling  between 
infragraviiy  waves  (ft  <0.04  Hz)  and  waves  with  frequencies  within  the 
power  spectral  peak  (0.04  Hz  +/2<0.12  Hz). 
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3).  in  shallow  {relative  to  a  sea  and  swell  wavelength)  water 
is  particularly  sensitive  to  changes  in  the  angular  spread  between 
interacting  free  waves  (Figure  2a).  In  O  (10  tn)  depth,  even  slight 
directional  spreading  of  the  sea  and  swell  spectrum  reduces  E^ 
to  10%  of  the  value  predicted  with  a  unidirectional  wave  field.  In 
contrast,  the  bound  wave  coupling  coefficient  in  deep  water  has 
the  largest  magnitude  for  large  angular  spread  (Figure  Ifi)  and  a 
directionally  broad  sea  and  swell  spectrum  can  force  more  bound 
wave  energy  than  a  relatively  narrow  one  (Figures  2c  and  2d). 
Directionally  spread  seas  in  deep  water  may  excite  "short"  (high 
wavenumber)  infragravity  waves.  Because  these  short  infragrav¬ 
ity  waves  will  not  penetrate  deeply  into  the  water  column, 
infragravity  energy  measured  at  the  sea  suifacc  in  deep  water  can 
be  radically  different  from  infragravity  energy  measured  on  the 
seafloor  or  even  a  few  hundred  meters  below  the  surface.  Bound 
wave  energy  is  also  a  strong  function  of  the  free  wave  frequen¬ 
cies.  In  shallow  water,  forced  by  low-frequency  (0(0.07 
Hz))  swell  is  10^  times  larger  than  forced  by  equally  ener¬ 
getic  high-frequency  (O  (0.2  Hz))  seas  (Figure  3).  In  deep  water, 
high-frequency  seas  force  more  E'*^  at  the  sea  surface  than  swell, 
although  the  differences  (sea  versus  swell)  in  E*^  energy  are 
smaller  than  in  shallow  water. 

£•*•  from  three  sites  in  8-  to  183-m  depth  were  compared  with 
predictions  based  on  unidirectional  and  directional  bound 
wave  theory.  Low-resolution  frequency-directional  incident  wave 
spectra  were  estimated  using  a  small  6  m  x  6  m  slope  array.  In 
-8-m  depth  at  Barbers  Point,  the  unidirectional  is  frequently 
greater  than  E’*^  and  exceeds  the  directional  E*^  by  as  much  as 
a  factor  of  20  (Figures  5,  lb,  and  Id).  Directional  spreading  of  the 
incident  wave  field  must  be  included  when  applying  bound  wave 
theory  in  shallow  water.  The  ratio  E'^/E’*’,  the  fraction  of  the 
infragravity  energy  which  is  bound,  varies  with  depth  and  sea  and 
swell  energy  {Ess).  The  portion  of  £“•"  accounted  for  using 
directional  bound  wave  theory  ranges  from  less  than  0.01  in  13-m 
depth  with  small  Ess  to  -0.4  in  8-m  depth  with  large  Ess  (Figure 
8).  In  183-m  depth  -10  km  offshore  of  Point  Conception,  uni¬ 
directional  is  less  than  10%  of  £‘*‘,  even  with  large  Ess 
(Figures  5  and  10).  The  estimates  of  £*^  are  qualitative  because 
of  uncertainties  in  the  sea  and  swell  directional  spread,  but  also 
because  in  shallow  depths  the  theoretical  assumptions  of  small 
Ursell  number  and  a  slowly  varying  wave  field  may  be  violated. 
However,  small  £’^/£“**  ratios  occur  in  all  depths  and  show  sys¬ 
tematic  variations  which  cannot  be  attributed  to  the  low  resolution 
of  the  sea  and  swell  directional  spectrum.  Free  infragravity 
waves,  either  leaky  waves  or  edge  waves,  are  more  energetic  than 
bound  waves  at  both  the  shallow  and  deep  sites.  The  low  level  of 
infragravity  energy  observed  in  183-m  depth  compared  with  8-  to 
13-m  depths,  with  similarly  moderate  sea  aixl  swell  energy,  sug¬ 
gests  that  leaky  (and  very  high  mode  edge)  waves  contribute  less 
than  10%  of  the  infragravity  energy  in  8-13  m.  Most  of  the  free 
infragravity  wave  energy  in  shallow  water  is  refractivciy  trapped 
and  does  not  reach  the  deep  site  (Figure  11).  At  all  of  the  sites 
and  depths  considered,  £^  and  £'*•  were  highly  correlated, 
suggesting  a  strongly  related  generation  mechanism  for  free  and 
bound  infragravity  waves.  Bound  waves  in  shallow  water  may 
indeed  be  a  source  of  free  infragravity  waves,  as  suggested  by 
Longuqt-Higgins  &  Stewart  [19^2]. 

1  “  '  ,  •  ■  ■ 
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The  first  two  Fourier  coefficients  a\{f),  b\{J)  of  the 
frequent-directional  spectrum  were  estimated  from  the  slope 
array  cross  spectra  [Ho'bers  and  Guza,  1989],  Following  Kuik  et 
al.  {198^1,  at  eacKftequency  the  mean  direction  6d(/ )  is 


0a(/  )  =  aiclan 


b^(f) 

_^K77J 


(Al) 


and  the  rms  angular  deviation  of  energy  from  the  mean  direction 
is  (for  narrow  distributions). 


cr(/)  = 


'f  (e-eofD(f  .e)de 
1  <  « 


=  l2(I-m!)r  (A2) 


where 


(A3) 


The  directional  distribution  was  modeled  as  the  widely  used 
parametric  form 


D{f  ,9)  =  C(/’)cos^</)| 


e-Odj) 


1 


(A4) 


where  C(f)  is  a  normalization  constant  and  the  spreading  param¬ 
eter  s(f)  was  selected  (Afuiil  el  al.,  1988,  equation  38]  so  that 
D{f  ,9)  (equation  (A4))  had  directional  width  equal  to  aif) 
(equation  (A2)): 


s(f)^ 


2 

W7 


(A5) 
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